Optimal decoding in fault-tolerant concatenated quantum error correction 



0^ 
O 
O 
(N 



Oh. 
a ■ 

^3 ; 



> 
o 

o 

ON 

o 



X 



Zachary W. E. Evans* and Ashley M. Stephens 
Centre for Quantum Computer Technology, School of Physics 
University of Melbourne, Victoria 3010, Australia 
(Dated: February 26, 2009) 

Here we present a new approach to concatenated quantum error correction in which additional 
classical processing is used with existing quantum codes and fault-tolerant circuits to more reliably 
correct errors. Using this new scheme we are able to guarantee correction of a number of discrete 
errors up to half of the distance of the concatenated code. We simulate the application of the scheme 
to the [[4,1,2]] subsystem code with up to three levels of error correction. 
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I. INTRODUCTION 

Our expectation that a quantum computer will one 
day outperform a classical computer is in part founded 
upon quantum error correction, which will be necessary 
to mitigate environmental decoherence and errors that 
arise from imprecise quantum control. Various methods 
for error correction are known, including topological er- 
ror correction [U, d, 0| , surface codes ^1, ^, ^] , and color 
codes @ . The method we are concerned with in this pa- 
per involves concatenating small quantum codes that are 
analogous to classical codes such as the repetition code 
and the Hamming code. 

Quantum codes such as the Bacon-Shor code and 
the Steane code [§] work by redundantly encoding the 
state of a single, physical qubit in the combined state of a 
few qubits, whereafter an error in the state of one or more 
of these qubits can be diagnosed and corrected using ap- 
propriately designed circuits that control the propagation 
of errors 10]. This reduces the error rate of an encoded 
quantum circuit relative to an unencoded quantum cir- 
cuit by at least a factor proportional to reciprocal of the 
physical error rate, as two or more concurrent errors are 
required to cause failure. However, because there are 
more places for errors to occur than before, the encoded 
error rate will only be lower than the physical error rate 
if the physical error rate is below some threshold [Tlj . 

A single level of error correction can only reduce the 
encoded error rate so much. A greater reduction can be 
achieved by encoding a single physical qubit in the state 
of a few qubits, as before, and then encoding each of those 
qubits in the state of more qubits, and so on. This forms 
a concatenated quantum code, which allows, in theory, 
an arbitrarily long and large universal quantum compu- 
tation to be undertaken reliably and efficiently. However, 
as both the number of qubits and the time required in- 
crease exponentially with the number of concatenation 
levels, it is likely that only a few levels will be practical. 
Therefore with finite resources the fault tolerant thresh- 
old should not be thought of as the critical point for 
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achieving arbitrary accuracy in quantum computing. To 
achieve a sufficiently low logical error rate with a practi- 
cal amount of resources it is important to make each level 
as effective as possible at reducing the encoded error rate. 

Typically, each level of a concatenated code operates 
independently of all other levels and implicitly assumes 
that errors on all qubits are equally likely. For a distance- 
d code that is concatenated I times, this method results 
in a concatenated code that can fail if there are |^rf/2]' 
errors. However, as the total distance of the concatenated 
code is c?', one might expect that a more sophisticated 
decoding method may exist which will never fail with 
fewer than [d'/^l errors. Such a method is known for the 
case of quantum information being communicated over a 
noisy channel [T^ . but not for the case where errors can 
occur during error correction itself. 

Recently we reported on a new method of concate- 
nated error correction which makes use of classical infor- 
mation obtained during the first level of error correction 
to assist in the diagnosis of errors at the second level 
p^ . When applied to the [[9,1,3]] subsystem code our 
method was shown to correct all cases with four physical 
errors whereas conventional error correction, where each 
level of error correction operates independently, can only 
guarantee success with up to three physical errors. Our 
method was inspired by a similar method that has been 
applied to an error-detection code with the result that 
the minimum number of errors that cause failure scales 
with increasing levels of concatenation as the Fibonacci 
sequence [13: [S] • 

Here we extend the method of Ref. [3] so that it is 
applicable to any number of levels of concatenation. By 
using information gathered by the process of error correc- 
tion itself we determine the relative likelihood of failure at 
each circuit location and thus maximize the probability 
of correctly interpreting high-level error syndromes. Our 
method requires no additional quantum resources and in- 
creases the minimum number of errors that cause failure 
to half the total distance of the concatenated code. In the 
design of this method we assume a discrete error model 
in which the probability of error at each physical circuit 
location, p, is constant and independent of the other lo- 
cations. We also assume that the quantum computer will 
be operating in the low-p limit, thus the exponent of the 
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lowest order of p is sufficient to describe probabilities. It 
may be possible to improve the method by relaxing these 
assumptions. 

II. METHOD 

The key idea for our method is that the process of 
error correction provides information that can be used 
to estimate the probability that an encoded qubit has 
failed. For example, whenever error correction measures 
a non-zero syndrome we know that at least one error 
has occurred, and so our confidence in the encoded in- 
formation should be decreased. Here we aim to use all 
of the information at our disposal to determine the most 
likely set of errors that could give rise to the measured 
syndrome so that we can maximize the probability of re- 
turning to the correct codeword state. We achieve this by 
tracking effect of each possible error through the error- 
correction circuitry, along with an associated confidence 
rating. Whenever a syndrome measurement is made, the 
results are compared to predicted syndrome measure- 
ments for possible sets of errors. We apply corrections 
to the data qubits based on the most likely set of errors 
that would give rise to the syndrome. Confidence ratings 
for the higher-level encoded qubits are calculated using 
the information from the syndrome measurements. 



A. Flag raising 

For each location in the error-correction circuit we de- 
termine what effect a discrete X and Z error occurring at 
that location would have on the data and ancilla at the 
end of syndrome extraction. Taking into account that 
some errors are equivalent to other errors, we assign a 
flag to each unique error. In other words, there is a flag 
for each possible effect on the data and ancilla that can 
be caused by a single error in the circuit. Each flag repre- 
sents the possibility of a single error. The probability of 
any given flag representing an actual error is described by 
the weight of the flag such that Pr(error) = ©(p'^^sht)^ 
As the possibility of error at each location is indepen- 
dent of the other locations, the probability of multiple 
flags representing actual errors is (^(p^ ^'shts-j^ where 
the sum is over the set of flags being considered. Error 
correction and flags in the X and Z bases are treated 
separately. 

At the physical level, each location in the error- 
correction circuit is assigned a confidence of one, indicat- 
ing that physical-level locations all fail with probability 
p. At the end of each error-correction cycle (in each ba- 
sis), after measuring the syndrome and applying the ap- 
propriate correction, a confidence rating is calculated for 
the encoded qubit based on the syndrome and flags. This 
is the logical confidence rating of the higher-level circuit 
location, the calculation of which will be described in the 
following section. Initially, all fiags can be thought of as 



having infinite confidence (as there cannot be any errors 
before any operations have been done). As each circuit 
location is executed, the weight of the flag corresponding 
to an error at that location is updated to be the minimum 
of its current value and the confidence reported at that 
location. Detection of errors at lower levels will lower 
our confidence of success at particular circuit locations 
all the way up to the highest level of error correction. 

Two-qubit gates are special in that errors may be 
shared on both qubits. For example, errors occurring 
before a transversal CNOT may be copied and errors oc- 
curring during a CNOT may affect both the control and 
the target qubits. Therefore, the confidence ratings re- 
ported by error correction immediately following a two- 
qubit gate are used to update the flag which corresponds 
to a correlated error affecting both qubits involved in the 
gate. This is in addition to updating both of the flags 
corresponding to the uncorrelated errors after gate. The 
flag weight corresponding to the two-qubit error is taken 
to be the maximum of the two single-qubit confidence 
ratings as opposed to the sum, which would correspond 
to the two uncorrelated errors. The possibility of a cor- 
related error is thus represented by its own fiag. 

Any error on a data qubit will eventually have its effect 
transfered down to the ancilla. The location of the error 
determines which ancillary qubit will see its effect. An 
error which affects a data qubit only after it has inter- 
acted with the ancilla in its current error-correction cycle 
can not be detected until the following error-correction 
cycle. Therefore, flags that represent errors that effect 
only data qubits should be put aside until the current 
error-correction cycle is complete, after which they will 
represent errors that will be detected in the following 
error-correction cycle. These flags can be called transi- 
tive flags. Following error correction, each transitive flag 
is converted into the its associated non-transitive flag for 
use in the next error-correction cycle. (See Section IIIII 
for an example.) 

Each logical qubit has its own set of flag weights. When 
two logical qubits interact (during a logical operation) 
certain errors within one logical qubit may propagate 
onto the other logical qubit. Let us call the logical qubit 
from which errors may be copied the source and the logi- 
cal qubit to which errors may be copied the target. Flags 
representing the possibility of errors on the source are 
used to update the weights of flags on the target. For 
each possible error that may be copied from the source, 
the corresponding flag weight on the target is taken to be 
the minimum of its current weight and the weight of the 
flag on the source. Although copying an error from one 
logical qubit to another would actually result in a corre- 
lated error affecting both qubits, we treat each error and 
its associated probability as if they were independent so 
that error correction can be performed on each logical 
qubit separately. The possibility of correlated errors is 
reintroduced at the logical level where the weight is taken 
to be the maximum of the weights of the individual log- 
ical qubits. 
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B. Flag matching, correction, confidence reporting 

At the end of each error-correction cycle the anciUary 
qubits are measured to obtain the syndrome. Each pos- 
sible syndrome is consistent with a finite number of sets 
of errors. Each of these sets of errors has a corresponding 
set of flags. We call each set that matches the measured 
syndrome a flag match. The most likely cause of the syn- 
drome is the set of errors that is represented by the flag 
match with the lowest weight, where the weight of a flag 
match is the sum of the weights of the flags in the set. 
With the lowest weight flag match identified, we apply 
corrections to whichever data qubits would be affected 
by the set of errors represented by the match. 

Alternative flag matches will imply different data cor- 
rections. In particular, there will always be a flag match 
which implies a complement set of corrections to that of 
the lowest weight match — that is, the corrections of both 
matches combine to form a logical operator. If we cor- 
rect based on some match but the true set of errors are 
represented by its complement match then the result is a 
logical failure. To determine the logical confidence rating 
that will be used by error correction at the level above, 
we consider the probability that our choice of correction 
will result in a logical error, 



Pr(fail|Sn flags) 



Pr(failnS|flags) 
Pr(S|flags) 



(1) 



Pr(fail|Snflags) is the probability of a logical failure given 
the current syndrome and set of flags. Pr(S| flags) is 
the probability of the current syndrome occurring given 
the current set of flags, regardless of success or fail- 
ure. Pr(failnS|flags) is the probability of the current 
syndrome occurring and error correction resulting in a 
logical error. In the low-p limit these probabilities are 
dominated by the leading-order terms and the coefficients 
are not important. Since our method applies corrections 
based on the most likely cause of the syndrome, to lead- 
ing order Pr(SS)=Pr(successnSS). This implies 



Pr(fannS|flags) 
Pr (success fl S| flags) 
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(2) 



where x is the weight of the lowest weight flag match and 
y is the weight of the complement match. Therefore, 



Cl - X 



is the logical confidence rating calculated during this 
error-correction cycle. 

The logical confidence rating represents the probabil- 
ity of a logical error, 0{p'^^). However, there is also a 
possibility that the corrections applied will result in a 
state that is outside of the code space, so that it is nei- 
ther the correct state nor the state affected by a logical 
error. Just as we used the complement match to deter- 
mine the probability of a logical error we can determine 
the probability of individual errors on each data qubit by 



considering the other flag matches. At the end of each 
error-correction cycle the confidence of each data qubit is 
updated to be the minimum of the confidence of its asso- 
ciated transitive flag and a confidence determined during 
the error correction of the encoded qubit, where this con- 
fidence is calculated in an analogous way to the logical 
confidence. 

In large codes there will also be probabilities of corre- 
lated sets of data errors that we are able to derive from 
the flag matches. Let E represent a particular set of a 
data errors after an error-correction cycle. Then 



Pr(S|Sn flags) 



Pr(i;nS|flags) 
Pr(S|fiags) 



(4) 



The weight of the flag corresponding to this set of data 
errors should, therefore, be updated to be the minimum 
of its current weight and the sum of each set of flags that 
match this outcome, taking into account any corrections 
that have been applied, minus the weight of the minimum 
weight flag match. (See Section Hill for an example.) 

By the careful consideration of every possible cause 
of every syndrome we always have accurate confidence 
ratings for every qubit at every level. We use this infor- 
mation to maximize the probability of success by always 
correcting based on the most likely set of errors. 



III. APPLICATION 

Here we apply our method to the [[4,1,2]] subsystem 
code 0, Q ■ The stabilizer generators of this code are 



X1X2X3X4, 
Z1Z2Z3Z4, 



(5) 



where Xi and Zi respectively represent the Pauli opera- 
tors (Tx and cTz applied to the i*'* qubit. Tensor products 
and identity operators are implicitly present. The gauge 
group of this code is generated by the operators 



XlX2, 

X3X4, 
Z1Z3, 
Z2Z4. 



(6) 



(3) The logical X and Z operators are 



Xl — A1A4, 
Zl = Z1Z4, 



(7) 



or equivalcntly any operators that can be generated by 
combining Eq. [7] with the gauge operators of Eq. [6l 

The syndrome-extraction circuit which we have chosen 
to use is shown in Fig. [1] [l^ , along with the flags that 
will be used in our message-passing method. The circuit 
performs operator measurements of the gauge operators, 
two at a time. Since the stabilizers are products of the 
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gauge operators, the combined parity of each pair of an- 
cilla measurements is effectively the result of an operator 
measurement of a stabilizer. Individually, each gauge op- 
erator measurement tells us nothing about the syndrome. 

Since the logical operators of the [[4,1,2]] code are pair- 
wise products of single-qubit operators, the code distance 
is equal to two. This means that it is only an error- 
detection code. The syndrome of the code consists of a 
single bit in each basis. A syndrome of zero indicates that 
no error was detected; one indicates that an error was de- 
tected but provides no information on which qubits were 
affected by the error. When the syndrome is measured 
to be one, in the absence of any additional information, 
we must simply guess which data qubit was affected by 
the error. Due to the gauge group of the code, each data 
qubit has a partner qubit in each basis for which errors 
are equivalent. For example, an X error on the first qubit 
is the same as an X error on the second qubit up to the 
application of the X1X2 gauge operator. Such an error 
can, therefore, be corrected by the application of an X 
operator on either the first or second qubit. To simplify 
the analysis of errors, we refer to errors by which gauge 
operator that they are a part of, rather than by which 
data qubit they affect. An X error on the first or second 
qubit is said to be an error on the first gauge in the X 
basis, a Z error on the second or fourth qubit is said to 
be an error on the second gauge in the Z basis, and so on. 
When there is an error on one gauge and the syndrome is 
misinterpreted so that the correction operation is applied 
to the other gauge, the result is always a logical error. 

In order to reliably correct errors using the [[4,1,2]] 
code, extra information is required to distinguish be- 
tween errors affecting the first and second gauges. One 
method for obtaining this information is to pass classi- 
cal messages describing the relative probabilities of er- 
rors from lower levels of error correction to higher levels. 
Previously, it was shown that a message passing scheme 
such as this could increase the minimum number of er- 
rors that cause failure to a number that scales according 
to the Fibonacci series with increasing levels of concate- 
nation [3, [IB]- In contrast, the scaHng achieved by our 
method is given by 2'~^ where I is the number of levels 
of error correction. 

It is worth noting that we have chosen to use [[4,1,2]] 
code for several reasons. Since it is a small code we 
can simulate more levels of error correction than if we 
had used a higher distance error-correction code. As 
our method involves an analysis of combinations of er- 
rors that can occur at every circuit location, it is easier 
to understand and test the operation of the method on 
a code that requires only small circuits with relatively 
few locations. Finally, a similar code forms the basis of 
the post-selection scheme of Knill, where ancillary states 
are prepared and verified before they are used to enact 
error correction by teleportation. The advantages of this 
scheme include a high threshold and fast logical gates, 
though unfortunately these come at a high resource cost. 
Our method might be used to improve the performance 



Flag 


Effects 


AGl 


first gauge and ancilla 


AG2 


second gauge and ancilla 


A 


ancilla only 


Gl 


first gauge only (transitive) 


G2 


second gauge only (transitive) 



TABLE I: Names assigned to the flags of the [[4,1,2]] code. 
There are flve flags for each basis. First gauge refers to an 
error on either of the two qubits involved in the flrst gauge 
operator of the appropriate basis in Eq.|S] similarly for second 
gauge. Ancilla refers to an error which flips the parity of the 
two gauge operator measurements — that is, an error which 
changes the outcome of the syndrome measurement. Flags 
Gl and G2 become AGl and AG2 respectively at the end of 
the error-correction cycle as they will be seen on the ancilla 
of the following error-correction cycle. 



and overhead of the scheme. 



A. Technique 

For the [[4,1,2]] code, using the circuits we have cho- 
sen, there are five unique effects that can result from a 
single error. Therefore there are five flags that we are re- 
quired to track for each logical qubit at each level of error 
correction. These flags are defined in Table I. Note that 
two of the flags are transitive. The set of error locations 
belonging to each flag are shown in Fig. [1] 

In our scheme the errors represented by the flag match 
with the lowest weight are identified as being the most 
likely cause of the observed syndrome. Corrections are 
applied to any data qubits marked by the flag match. 
For example, if the X basis syndrome is one, and the 
weights of the X flags are AGl = 4, AG2 = 0, ^ = 2, 
then the flag match with the lowest weight is AG2 and 
so a correction should be applied to the second gauge — 
that is, an X correction should be applied to either qubit 
three or qubit four. Table II shows a list of all possible 
flag matches along with the various confldence ratings 
that would result from each of them being acted on. 



B. Simulation and results 

Our simulations are of a CNOT extended rectangle with 
concatenated error correction, including all physical lo- 
cations. The message-passing scheme is simulated along 
side the concatenated error correction circuit so that it 
operates in the same way as it would in a real quan- 
tum computer. The logical CNOT extended rectangle is 
defined to succeed only if both logical qubits would, af- 
ter decoding, produce the correct measurement results in 
both the X and Z bases. The simulations use the method 
outlined in p^ . which involves keeping track of discrete 
Pauli errors as they propagate through a circuit. We 
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AGl 



AGl 



AG2 



AG2 



1+) 
1+) 



N 



N 



AGl 



N 



AG2 



N 



X 



Gl 



A 



X 



Gl 



G2 



G2 



FIG. 1: Syndrome-extraction circuit for the [[4,1,2]] code, with X flag mapping shown. The error syndromes in the X and 
Z bases are obtained by performing operator measurements of the gauge operators in Eq. |6] The parity of the two X gauge 
operator measurements gives the Z error syndrome and vice versa. The letters on the circuit mark out regions associated with 
each of the X flags. The confidence rating of each circuit location is associated with the flag that is labeled to the left of that 
location. (See Tab. I for the flag labels). For example, an X error that occurs in the region marked ^Gl will affect one of the 
qubits involved the first X gauge operator and the effect will be seen on the ancilla used to perform the operator meeisurement. 
An X error that occurs in the region marked A will fiip the parity of the ancilla but will not affect the data at all. An X error 
that occurs in the region will have no effect on the data or the observed syndrome. The flag regions shown here are for X 
flags only but the Z flag map is similar. 



Syndrome 


Flag match 


Correction 


Gai 


Gg2 


Gl 




None 


None 


AGl + A 


AG2 + A 


AGl + AG2 





AGl + A 
AG2 + A 
AGl + AG2 




N/A 
N/A 
N/A 






AGl 


Gl 


A- AGl 


AG2 + A 


AG2 - AGl 


1 


AG2 


G2 


AGl + A 


A-AG2 


AGl ~ AG2 




A 


None 


AGl- A 


AG2-A 


AGl + AG2 




AGl + AG2 + A 




N/A 





TABLE II: Table of all flag matches for each of the two possible syndromes. Corrections and flag updates are based on the flag 
match with the lowest weight from the list in the second column. The correction and flag update rules associated with each 
flag match are shown in the last four columns. Gai and Gg2 are confldence ratings for the transitive flags Gl and G2 at the 
current level of error correction. Gl is the logical confidence rating for the next-highest level of error correction. 



assume a stochastic error model, where memory, initial- 
ization, readout, single-qubit gate, and two-qubit gate 
errors are all equally likely to occur. Errors are ran- 
domly selected one- and two-qubit discrete Pauli errors. 
Our simulator uses the SIMD-oricnted Mersenne Twister 
pseudo random number generator [17| . We use the error- 
correction circuit shown in Fig. [1] which requires two 
ancilla qubits per logical qubit for a total block size of 
six qubits. 

We are primarily interested in the logical failure rate 
when the physical error rate, p, is low. Directly simulat- 
ing the circuit in the low-p limit requires too many runs 
to obtain statistically significant results. Instead, we use 
the expansion of logical failure rate with respect to the 
physical error rate, 

PL-T.rJ^)p\l-pf-\ (8) 

where A'' is the number of physical locations in the error- 
correction circuit, p is the physical error rate, and is the 



probability of logical failure after error correction given 
exactly i errors. To estimate values of we simulate the 
circuit with exactly i errors placed randomly, repeating 
many times for each i. We approximate pl by truncating 
the series above i — 29. This approximation breaks down 
for values of p for which the probability of having more 
than 29 errors becomes significant. 

The failure rate of the logical CNOT as a function of 
physical error rate is plotted in Fig. [5] for one, two, and 
three levels of error-correction. In the low-p region, the 
gradient of each line on the log-log plot signifies the min- 
imum number of errors that are required for the circuit 
to fail: one, two, and four errors for one, two, and three 
levels of error correction respectively. 

IV. FURTHER WORK 

As described, our scheme assumes that each physical 
location is equally likely to fail. However, this need not 
be the case. Relative probabilities of failure for Individ- 



6 




no FC 
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2 levels 
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10"^ 10"^ 10"^ 10"'* 10"^ 



FIG. 2: Failure rate of CNOT extended rectangle as a function 
of physical error rate for one, two, and three levels of er- 
ror correction with the [[4,1,2]] code and our message-passing 
method. 2a of statistical error above and below the mean is 
shown for each set of data. The minimum number of discrete 
errors that can cause failure is one, two, and four for one, two, 
and three levels of error correction respectively. This can be 
inferred by the gradients of the lines in the low-p region. The 
asymptotic threshold appears to be around 2 x 10~*. Since 
the first level of the code is unable to reliably correct any er- 
rors, the failure rate for one level of error correction is greater 
than the physical failure rate. 

ual physical locations are naturally represented in our 
scheme by the confidences at the physical level. That is, 
rather than setting the confidence of all physical loca- 
tions to one, we can choose to assign difi^erent confidence 
ratings to each type of operation or each individual phys- 
ical device. These confidence ratings could be based on 
data obtained by characterization of the gates. Alter- 
natively, data collected during error correction could be 



fed back to dynamically update physical confidence rat- 
ings in a similar way to how our method updates logical 
confidence ratings. 

Our message-passing method could be applied to a 
post-selection based scheme instead of the Fibonacci 
scheme. There are several ways this could be done, since 
we are free to choose the confidence threshold which 
determines whether ancillary states are accepted or re- 
jected. To only accept states that have the maximum 
confidence would reduce the scheme to ordinary post- 
selection, as a state would be rejected if any errors are de- 
tected at any level regardless of whether or not these er- 
rors are correctable at higher levels. Accepting all states 
that have a confidence greater than one would result in 
a high acceptance rate, but would not be optimal with 
respect to logical fidelity. In any case, the confidence of 
the states that are accepted could be used in higher-level 
error correction following post selection. There may also 
be some advantage to using the logical confidence ratings 
outputted by the highest level of error correction in some 
quantum algorithms. 

Our method can be applied to higher distance Bacon- 
Shor codes and perhaps even to other CSS stabilizer 
codes. It would be interesting to investigate the pros and 
cons of using different codes with our method and also 
to rigorously establish the specific requirements on codes 
and circuits that enable us to achieve optimality with re- 
spect to the minimum number of errors that cause the 
code to fail. 
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